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Abstract 

Under periodic boundary condition in the transverse direction, we calculate 

the averaged zero-temperature two-terminal conductance (< G >) and its 

statistical fluctuations (< > for n < 4) at the critical point of integer 

quantum Hall plateau transitions. We find universal values for < G >= 

(0.58 ±0.03)4, < i^G?" >= (x)^"^2n, where ^2,4,6,8 = 0.081 ± 0.005; 

0.013±0.003; 0.0026±0.005; and (8±2) x lO""^ respectively We also determine 

the leading finite size scaling corrections to these observables. Comparisons 

with experiments will be made. 
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For macroscopic, disordered, two dimensional electronic systems at low temperatures, 
the metallic behavior is only observed near quantum phase transitions. Two examples are 
the superconductor to insulator transitions in amorphous thin films, and the transitions 
between quantum Hall plateaus In two spatial dimensions, the conductivity tensor 0"^,^ 
measured in units of e^/h is dimensionless. Under generic conditions, < cr^j, > (the impurity 
averaged conductivities) are expected to be universal at the quantum critical points 

The universal conductance fluctuations of disordered mesoscopic metals has attracted 
tremendous experimental and theoretical interests in recent years 0]. The same question 
can be asked at the quantum critical points mentioned above. Thus "what is the statistical 
(over the impurity ensemble) properties of a^i, at quantum critical points of two dimensional 
systems" is the concern of the present paper. In particular, we concentrate on the transitions 
between integer quantum Hall plateaus. 

In an integer plateau transition, the electron conductivity tensor [axx, o'xy) changes from 
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(0,n)^ to {0,n ± 1)^. Concerning these transitions the following consensus are reached 
IHIJl^- (^) When extrapolated to zero temperature and inflnite sample size, the transitions 
are genuine continuous phase transitions with a single divergent length scale ^ - the quasi- 
particle localization length, (b) As the Fermi energy (Ep) moves toward a critical value 
Ec, ^ ~ \Ef — Ec\^'^, with u ~ 2.3. (c) The characteristic energy scale ~ 1/^, thus (the 
dynamical exponent) z = 1. 

Based on a Chern-Simons formalism, Kivelson, Lee, and Zhang (KLZ) asserted that 
i<Jxx,<Jxy) at the (0,n)^ (0,n + 1)^ transition is (a^x^^xy) = (|?^ + |)x 0- 
Hetzel, and Bhatt have numerically computed a^^ at the (0,0) (0,1)^ transition using 
the Kubo formula 0. Their calculation assumes that the electronic states lie entirely in the 
lowest Landau level. They have considered several different forms of the impurity potential 
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and concluded that a^^ = (0.55 ± 0.05)^ which is consistent with the result of KLZ. 

Recently, the experimental four-terminal resistivity pxx has been studied systematically 
at the quantum Hall liquid to insulator transitions by Shahar et. al. [^]. They concluded 
that the critical conductivity tensor is universal and the value of p^^ h/e^ agrees with 
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the prediction KLZ. In addition, significant sample to sample fluctuations in the critical 
resistance have been observed [0]. Here we note that in order to measure the true p^^ 
(or cr^a;), the temperature has to be low enough so that corrections to scaling have faded 



away, yet it has to be high enough so as to avoid the effects of mesoscopic fluctuations [11 



More recently, the statistical fluctuation of the two-terminal conductance in the transition 
regime has been studied experimentally in mesoscopic samples by Cobden and Kogan who 



demonstrated the presence of large mesoscopic fluctuations in the conductance [12 



In this paper, we calculate the ensemble averaged two-terminal conductance < G > and, 
for the flrst time, its fluctuations < > at the critical point of the integer quantum 

Hall plateau transitions. In addition, we study the flnite-size and aspect ratio dependence 
of these quantities. The model we use is the Chalker-Coddington network model and its 



extension [|1^], and with periodic boundary conditions applied in the transverse direction. 
These models have been shown to exhibit the correct critical properties of the integer plateau 
transitions |r^-p!5|. Our flndings are summarized as follows. The conductance of a x L 
sample (where W is the circumference and L is the length of the cylinder) exhibits the 
following scaling behavior, 

<G>= jT{Ly-^Ag,,u L-y-^Ag,,,, W/L). (1) 

Here we have used the fact that G/{^) is dimensionless. In the above, < ... > denotes 
the disorder ensemble averaging, Ag^^i and Ag^j-j. are the coupling constants conjugate to 
the relevant and the least-irrelevant operators respectively. Thus i/rei = 1/z^ and — are 
the renormalization group dimensions of Ag^-ei and Agi„ respectively. At the critical point, 
= 0. We have performed finite size scaling analysis of< G >c on L x L samples 
to extract JF(0,0, 1) and to be (0.58 ± 0.03) and 0.55 ± 0.15 respectively. Thus the 
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critical conductance < G >c= (0.58 ± 0.03)^. We have also calculated the central moments 
< (5G)^" > for n < 4 and shown that they exhibit the following scaling behavior, 

< 6G^^ >= (^)2"^2n(i^^-'A^,ei, L-y^-Ag,,,, W/L). (2) 
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We have determined the critical moments < 5G'^^ >c= (e^//i)^"'JF2„(0, 0, 1). For n = 1,2, 3, 4, 
the values of J^2n are found to be 0.081±0.005, 0.013±0.003, 0.0026±0.005, and (8±2)xl0-^ 
respectively. We verified the universality of these results and assert that all higher moments, 
and hence the entire distribution function P{G), is universal at the transition. Since in the 
rest of this paper we ignore the electron-electron interaction, our conclusions are at most 
relevant to the integer plateau transitions. We also point out that it remains unclear at 
present about the relation between the two-terminal conductance G and the a^^,^ derived 
from the Kubo formula in a closed system without contacts. 

Let us consider non-interacting electrons in a strong magnetic field and potentials that 
are smooth on the scale of the magnetic length. In this limit the plateau transition can be de- 
scribed by "quantum percolation" of semiclassical electron orbits in the Chalker-Coddington 
network model |[13|JT^ . The latter consists of a square lattice of potential saddle points at 
which quantum tunneling between the edges of quantum Hall droplets take place (Fig. la). 
Away from these vertices, the edge electrons propagate along the directed links with a 
fixed chirality set by the direction of the magnetic field. To account for the random areas 
of the Hall droplets, the edge electrons accumulate random Bohm-Aharonov phases while 
traversing the links of the network. At each saddle point, as shown in Fig. lb, there are 
two incoming and two out-going edges states. The associated probability amplitudes are 
denoted by ^i,...,4. Due to current conservation, + 1^2^ = 1^3^ + |^4p. With a choice 
of gauge, the quantum tunneling event at each node is then completely specified by a 2 x 2 
transfer matrix 

ZA fe'^^ \ /cosh^ sinh^\ ( e'^^ \ f 
Zs) V e^<^2 j y sinh 6 cosh ^ / V e'^^ / V ^2 
where (pi G [0, 27r) are random phases, and 6* is a real number. Using Eq. as the building 
block, we construct the W xW transfer matrix Tj which propagates eigen wavefunctions on 
a cylinder of circumference W in the x-direction. The total transfer matrix T for the entire 
W X L system is given by the matrix product T = Hi^i Ti- For the details of transfer matrix 
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calculations readers are referred to the original work of Chalker and Coddington |]13[. Here 



we merely emphasize a few important points. 

In general the tunneling parameter 9 in Eq. should be a random variable depending on 
the local potential at the saddle point. However, in Ref. [jl4[ it was shown that introducing 



randomness in 9 did not change the universality class [i.e. did not change the localization 
length exponent v = 2.33 ± 0.03) of the plateau transition. It turns out that the network 



model is situated at the critical point as long as 9 is the same for all nodes [1^. Invariance 



under 90-degree rotation selects a particular one (9 = 9c = ln(l + v^)) out of the family 



of critical models ||I3|,|15|. We shall present results for both isotropic systems {9 = 9c) for 
which the two-terminal conductance along the x-direction (Gxx) and the y-direction {Gyy) 
are equal, and anisotropic systems (i.e. 9 ^ 9c) for which G^x 7^ Gyy. In the latter case we 



study the geometric mean ^jG^xGyy. 

We next define the two-terminal conductance. As shown in Fig. la, two semi-infinite 
conducting leads are connected to the W x L disordered network. It has been shown that 
linear response theory applied to the combined lead-sample-lead system gives the following 
multi-channel two-terminal conductance [ffB,0], 



p2 2 W 1 

G = ^Trttt = ^V i -. (4) 

h hfr[ cosh2(7iL) ^ ' 

In Eq. (^), t is the transmission matrix through the disordered region and 7j is the z-th of 
the W Lyapunov exponents of the hermitian transfer matrix product T'^T [|l^]. We note 
that the validity of Eq. in the absence of time reversal symmetry has been shown in 
detail by Baranger and Stone for finite magnetic fields 

We now present the results. We have studied W x L systems with periodic boundary 
conditions in the transverse direction for L = 8, 16, 32, 48, 64, and 96, and for aspect ratio 
W/L = 1/4,1/3,1/2,1,2,4. The disorder ensemble consists of at least 3000 samples for 
each {W,L). The distribution function P{G) is very broad, and the most probable value 
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of the critical conductance Gctypicai is, although close to 0.5^, not sharply defined. For 
isotropic systems (i.e. 9 = 9c) the averaged critical conductance < G{L) >c for L x L 
samples are shown in Fig. 2. To extract the asymptotic value < G >c and the exponent 



?/irr at the critical point, we expand the scahng function in Eq. (|l|) according to JF(0,x, 1) ~ 
^-"(0, 0, 1) + J-''{0, 0, l)a; for small x. Thus for large system size L we should have 

< G{L) >e=< G >, + J-'(0, 0, l)Ag;,,L-y-% (5) 

where < G >c= J-'{0,0,l)e'^/h. From the data shown in Fig. 2, we obtain JF(0,0, 1) = 
0.58 ±0.03 and i/ij-r = 0.55 ±0.15. We have also studied the aspect-ratio (W/L) dependence 
of < G >c- In the context of Eq. (|l]) we find 

^(0, 0, W/L) = cie-^2^/^(iy/L), (6) 

with Ci = 0.72 ± 0.03 and C2 = 0.22 ± 0.02. The exponential factor is a precursor to the 
conductance behavior in a quasi-one dimensional system with L ^ W. 

In order to verify the universality of < G >c and to check that the — ?/irr is indeed the 
renormalization group dimension of the least-irrelevant operator, we have studied network 
models with various forms of node parameter disorder corresponding to random distributions 
of ^-values with the same median 6c- The results are consistent within the error bars with 
those obtained above. This further supports the conclusion that ^-disorder is an irrelevant 
perturbation at the critical point. In addition, we have also considered anisotropic critical 
models for which 6 ^ 9c- In this case, using the procedures described above, we have 
calculated the conductances < G^x > and < Gyy > (periodic boundary conditions are 
always imposed in the transverse directions). 

The results are summarized in Table I. Although the critical values < G^x >c and 
< Gyy >c depend on the amount of anisotropy, their geometric mean ^GxxGyy does not. The 
latter stays close to the value of < G >c obtained for the isotropic system. Thus we conclude 
that there exists a universal critical two-terminal conductance < G >c= (0.58 ± 0.03)e^//i. 

The fact that our value for < G >c deviates from a^^ = asserted by KLZ is, as far 
as we can see, real. However, as was pointed out at the beginning of this paper, it is not 
clear how the two-terminal conductance relates to the conductivity tensor derived from the 
Kubo formula in systems without contacts. 
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We now turn to the central moments of the critical conductance described by the scaling 
forms given in Eq. (Q). In the inset of Fig. 2 we plot Ti =< >c /{e^/hY for the 

isotropic network (i.e. 9 = 6c) as a function of the system size for L x L samples. Following 
a similar finite size scaling analysis as given in Eq. (|^), we extract JF2(0,0, 1) = 0.081 ± 
0.005 and the same yirr = 0.55 ± 0.15 to be the renormalization group dimension of the 
leading irrelevant operator. The dependence of < (^G)^ >c on the degree of anisotropy 
{i.e. 9 ^ 9c) is shown in Table I. Thus we conclude < (SG)^ >c= (0.081 ± 0.005) (^)2. 
Repeating the procedure for the 4-th and the 6-th moments give the results < (5G)^ >c= 
(0.013±0.003)(^)^ < {5G)^ >c= (0.0026±0.005)(^)^ and with less accuracy < {5G)^ >c= 
(8 ± 2) X 10~^(^)'^. In Fig. 3, we present the data obtained for even-integer as well as odd- 
integer and non-integer central moments. The n-th order moments interpolate between 
< (^G)" >c= af"e""^ at small n and < (5G)" >c= hn~^ at large n, where {a,v,u,b, P) ~ 
(0.80,0.28,0.054,4.65,4.18) are all universal constants. The small-n behavior shows that 
the critical conductance obeys a log-normal distribution. On the other hand, the large- 
n behavior, indicative of a broad distribution, is purely empirical. The latter should be 
contrasted to the behavior in mesoscopic disordered metals where the the 2 + e-expansion 



in the diffusive regime predicts nonuniversal moments for large n ||T^ . 

We now compare our results with experimental findings. This comparison must be made 
under the disclaimer that so far the understanding of the effects of Coulomb interaction 



on the critical properties of the integer plateau transition is still incomplete pO[. First, 



we assume that the two-terminal conductance < G >c is the four-probe a^^ measured 
experimentally. Second, from general particle-hole symmetry argument one expects a^y = 
0.5^. With our result < G >c~ 0.58^, this implies = 0.99^. This value agrees with the 



experimental finding of Shahar et. al. |Ty] . More recently, Cobden and Kogan have measured 



the two-terminal conductance near the integer transitions. . Since the latter was carried 
out on mesoscopic samples outside the asymptotic scaling regime, direct comparison with 
our results would be difficult to justify. Nevertheless, their data provide an estimated value 
< (5G)^ >~ .05{e^ / h)"^ , which is in reasonable agreement with our findings. 
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FIGURES 




Figure 2. The finite-size dependence of tlie critical conductance at tlie isotropic critical 
point. The solid line is the fit to the scaling form in Eq. in the text with < G >c= 
(0.58 ± Om)e^/h and Vivv — 0.55 ih 0.15. Inset shows the behavior of the conductance 
fluctuations, which gives < {50^ >c= (0.081 ± 0.005) (eV/i)^ and the same renormalization 
group dimension i/i^^. 
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Figure 3. The behavior of the n-th order central moments of the conductance distribution 
which interpolates between the exponential and the power-law behaviors (dashed-lines) in 
the small and large n limits as described in the text. 
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TABLES 

TABLE I. Ensemble-averaged conductances and conductance fluctuations at the critical points 

of the isotropic {9 = 6^) and the anisotropic {6 ^ 9c) models. The units are /h for G and (e^//i)^ 
for {SGY . Error estimates in the last digits are given in parenthesis. 
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